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^ \ Abstract 

Solutions of gravitational equations of gauge theories of gravity in homogeneous isotropic 
world with massive scalar field are investigated in the case of flat cosmological models. Special 
attention is dedicated to general behavior of solutions on contraction stage. It is shown, that 
on expansion stage inflationary solutions are generic feature of the model. At the same time 
on contraction stage an effective equation of state is similar to the case of massless scalar field 
p = p at the beginning of evolution and tends to the equation of ultrarelativistic gas p = p/3. 
The Hubble parameter tends to some negative value on contraction stage, depending on the 
mass of the scalar field. Nonsingular solutions in this model are unstable. 
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Inflationary cosmology PQ-[E] allows to resolve a number of puzzles of Big Bang model. It became 
a standard paradigm since observational data favor a flat cosmological model and a scale-invariant 
| spectrum of initial perturbations, that are main predictions of inflation [S]. 

Cosmological singularity is a natural attribute of inflationary models in General Relativity 
o5 ' (GR). In spite of the fact, that nonsingular inflationary solutions exist in the case of closed models 
|0],[Zj, they cannot be used in order to construct viable nonsingular cosmological models. Their 
measure is equal to zero, and, moreover, on the contraction stage such solutions are unstable and 
have the effective equation of state p — p 

When departure from a Riemannian spacetime is allowed, the geometry becomes more general, 
such as afflne geometry with torsion and nonmetricity in addition to curvature tensor. A review 
of gauge theories of gravitation (GTG) in such spacetimes can be found in [E]. This could happen 
for example during the inflationary stage where in addition to the metric tensor a scalar dilaton 
induced by a Weyl geometry appears. 

One of the most important feature of GTG is a possibility to avoid cosmological singularity 
|lUj . Indeed, nonsingular cosmological models were built on the basis of GTG jTTj. In particular, 
inflationary nonsingular cosmological solutions were obtained ,12 : . 

However, the question on the generality of inflationary and (or) nonsingular solutions were not 
adequately explored within the framework of GTG in comparison with the case of GR (see for 
a review). In this paper solutions of cosmological equations of GTG in the case of flat models in 



1 



the Universe, filled by massive scalar field are analyzed and discussed. General behavior of such 
solutions is investigated both on contraction and expansion stages by means of qualitative theory 
of dynamical systems. 

Cosmological equations of GTG contain indefinite parameter (3 (see below) , that has a dimen- 
sion of inverse energy density and determines roughly the scale, when predictions of GTG diverge 
from that of GR. Notice, that regular cosmological models were obtained in the case (3 < 0. 

The paper is organized as follows. At the section 2 equations for a dynamical system are 
derived. In section 3 qualitative analysis is carried out. In section 4 the phase space at infinity is 
analyzed. Section 5 contains discussion. Conclusions follows in the last section. 



2 Basic equations 

Homogeneous isotropic models in GTG are described cosmological equations ^0] 
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where R(t) is the scale factor in the Robertson- Walker metric, k = —1,0,-1-1 for open, flat and 
closed models respectively, M p is planckian mass, p(t) and p(t) are energy density and pressure, 
f3 is indefinite parameter, and a dot denotes the differentiation with respect to time. 1 

Besides equations (0121) gravitational equations of GTG lead to the following relation for torsion 
function 5* and nonmetricity function Q 



1 d , 
— —ml 

Adt 1 



(3(p-3p)}. 



(3) 



Note, that in Poincare gauge theory of gravity Q = and equation (j^J determines the torsion 
function. In metric-affme theory of gravity equation (j3J) may describe three kinds of models: in the 
Riemann-Cartan spacetime Q = 0, in the Weyl spacetime S — 0, in the Weyl-Cartan spacetime 
(S 7^ 0, Q 7^ 0), the function S is proportional to the function Q 

The conservation law has the usual form 



p + 3H(p + p) = 0, 
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where H(t) is a Hubble parameter. 

In the case of massive scalar field minimally coupled to gravity the energy density and pressure 
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where m is the mass of the scalar field (p(t). 

hereinafter the system of units with h = c = 1 is used. 
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Then taking into account (jUEJ) the cosmological equations (012)) can be rewritten in the fol- 
lowing way 
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where Z = 1 + /3 [0 2 - 2m 2 p> 2 }. 

The conservation law (j3J) in the case under consideration is reduced to the scalar field equation 

<p + 3H<p + m 2 <p = 0. (8) 

In the case of flat cosmological models (k = 0) the Hubble parameter H in (JEJ) can be expressed 
in terms of {(p, <p}. Therefore, in the case k = the 3-dimensional dynamical system can be 
reduced to 2-dimensional one, by the analogy with GR Indeed, since H depends only on 
scalar field and its derivative, we have the following dynamical system: 
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where ip = tp. 

However, the expression for if in GTG contains a square root, that could become negative 
unlike the case of GR (the corresponding square root is always nonnegative). Consequently, phase 
space is divided into two parts, namely physical region, when scalar field is real and non-physical 
regions, when scalar field contains imaginary part. Condition 



/3(8m b /V + V~ PV) + W ~ 4m 4 /?</(2 + 3f3ip 2 ) + 2mV(l + M + 3(3 2 ^) = (10) 

determines the edge of the physical regions of the phase space. It is a sixth order algebraic 
equation and it provides the following three quadratic relations between dynamical variables: 
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ip 2 = /T 1 + 2/3mV ± ^1 + 6/3mV 2 ) • (12) 

From now on, we assume (3 < 0, since in the case of scalar field nonsingular solutions were 
obtained only for negative (3 ■ I n the case j3 < the edges of the physical region are determined 
by fTTf only. They are hyperbolas (see figH} with asymptotes 

ip = ±V2mip. (13) 
Notice, that asymptotes do not depend on the value of parameter f3. 




Figure 1: The edges of the physical phase space. Above and below upper and lower hyperbolas 
correspondingly dynamical variables contain imaginary parts. The tilt of asymptotes does not 
depend on the value of (3 and depend only on m (see (JUJ)). 

It is possible to show, that in the case (3 < sign "+" in (jHJ) corresponds to positive Hubble 
parameter, and, sign "-" takes place for negative Hubble parameter. Therefore for analysis of 
stage of expansion we will take equations Q with sign " +" , and for contraction stage we will use 
sign "-" in ©. 

3 Qualititive analysis 

We will formally consider infinite phase space, i.e. phase trajectories could go to infinity. Notice, 
however, that from the physical point of view, when the energy density (jSJ) becomes larger than 
the planckian value (p > M*), quantum gravity is needed. Therefore clear physical interpretation 
is possible only for those trajectories, that lay inside the quantum edge p = M^. 

In the finite region of the phase plane {<p, ip} there is only one singular point, namely the 
point {0, 0}. In the limit {(/? — > 0, ip — > 0} the system © is reduced to the case of GR. Therefore, 
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all consequences of the analysis of the corresponding dynamical system, that is given in [7j, take 
place in the system under consideration. 

Phase trajectories are shown at figlU Quantum edge (dashed line) marks the distance, at 
which the energy density is equal to the planckian value. Beyond the quantum edge classical 
interpretation of solutions becomes impossible. The edge of the physical region could come inside 
the quantum edge if \/3\ > M~ 4 (see (fTTjO. 
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Figure 2: Phase trajectories of flat cosmological models with massive scalar field within framework 
of GTG. Dashed line marks the quantum edge (see text). The picture corresponds to expansion 
stage. 

On the stage of expansion trajectories, starting at some particular point {<p , ipo} tend to 
inflationary separatrix, which lays close to the axis ip. Since at the inflationary separatrix ip <C <p 
(and m <C M p ) the effective equation of state is p ~ — p, i.e. inflationary stage takes place [Zj. 
Therefore, most solutions with particular initial conditions are inflationary. After inflationary 
stage the field oscillates near the point {0, 0} with damped amplitude. All phase trajectories 
finish at the spiral equilibrium point {0, 0}. 

The difference between GR and GTG cases is the presence of the edge of the physical region. 
Initial conditions cannot lay beyond these edges. On the stage of expansion they are repulsive 
separatrices. 

On the opposite, on the stage of contraction the inflationary separatrix is repulsive, and most 
initial conditions lead to the effective equation of state p ~ p. It is well known, that within the 
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framework of GR cosmological models with such equation of state have singularity in future at the 
stage of contraction. However, since in the GTG case there are two edges of the physical regions, 
the subsequent evolution of the model changes significantly. 

The edges of the physical region become attractive separatrices at the stage of contraction. 
Most trajectories finish at these edges. Notice, that since the edges are hyperbolas, i.e. infinite 
curves, dynamical variables tend to infinity in the course of time. 

Since the main difference between GR and GTG cases occurs at the contraction stage, below 
we consider H < everywhere. Note, that at the stage of contraction fig|2]must be reflected with 
respect to the vertical axis. 



4 Phase trajectories at infinity 

The asymptotic behavior can be characterized by the effective equation of state. Taking into 
account (J3J) and (fT3j) one arrives to 
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Thus, within the framework of GTG the effective equation of state for cosmological models 
with scalar field is the equation of state for ultrarelativistic gas 2 . Note, that in the case of equation 
of state (fllj) the cosmological equations (jl|2|) reduce to Einstein equations of GR. 

Hubble parameter can be calculated when the scalar field reaches the separatrix The 
corresponding value is 
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Therefore the Hubble parameter becomes constant and relatively small, taking into account 
that typical values of the mass are m ~ 10~ 6 M P jHj. In infinite future the contraction with constant 
speed will lead to singularity. Most trajectories finish at this state on the stage of contraction. 

To find all singular points of the system © we use the standard technique of projection of 
infinite space into the lower hemisphere of the sphere with unit radius, so called Poincare sphere. 
First of all, we transform cartesian coordinates {ip, ip} into polar, then project polar coordinates 
into the hemisphere and, finally, transform back to cartesian coordinates. 

The resultant transformations are 



u 



V 



■p 



arctg (i/j 2 + (p 2 )'- 



arctg [ip 2 + <p> 2 ) 2 



(17) 



2 The same equation of state appears also for nonlinear scalar field with effective potential V = jXip 4 and we 
believe that it is quite a common feature. 
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Figure 3: Phase portraits of the flat cosmological models with massive scalar fields in the cases 
of GR (left figure, see 0) and GTG (right figure), projected on the Poincare sphere. The dotted 
circle corresponds to the infinite values of the dynamical variables. The contraction stage is shown, 
H < 0. In order to get phase portrait at the stage of expansion figures have to be reflected with 
respect to the vertical axis. 

The result 3 is illustrated at figOl where left figure corresponds to the GR case, and the right 
figure corresponds to the case of GTG. The case of GR was investigated at paper [Zj and here we 
provide the corresponding figure for comparison. 

There are two knots, two saddles and one spiral focus (that we discussed in the previous 
section) in both cases. In the GR case knots lay at the axis u = and saddles lay at the axis 
v = at infinity of the real phase space. 

In the case of GTG two knots lay at the point, where asymptotes (fT3~j) cross the radius of the 
circle. Notice again, that the tilt of asymptotes does not depend on parameter (3. Moreover, in 
new variables the form of equation (|13|) is conserved. 

It could be seen, that the tilt of asymptotes ()13|) tends to zero with zero mass, i.e. it is possible 
to make the hyperbolas as flat as we wish. Taking into account that at the stage of contraction 
the hyperbolas ([lip are attractive separatrices, the question arise, what is the effective equation 
of state of the scalar field on these hyperbolas. 

From (0) and (fTTj) we have 

y = P (l + 3(l-3faV) )' (18) 

Therefore 

3 In effect, the radius of the circle, that corresponds to transformations i|17|) is it/2. 
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^<P<P- (19) 

The two limiting cases take place when the scalar field and its derivative with respect to time 
tend to infinity and to zero correspondingly. 

5 Discussion 

It was argued in [Hi], [IE] that it is possible to built nonsingular inflationary cosmological models 
within the framework of GR. However, in spite of the fact that such models could exist, according 
to the analysis of [7j their measure is equal to zero, because on contraction stage effective equation 
of state tends to p = p, that lead to cosmological singularity for most solutions. 

According to ^B], the possibility to avoid cosmological singularity remains if "physical condi- 
tions leading to the bounce take place at the end of cosmological compression" . It should be clear, 
however, that in this case nonsingular solutions have much less attractive character, because they 
necessarily require special physical conditions, for instance phase transitions, at the bounce. Such 
models are already discussed in the literature [T7|.[T8|. 

The analysis presented in the paper shows that the similar situation takes place within the 
framework of GTG in the case of flat models. Since cosmological equations of GTG are much 
more complicated than than of GR, they possess a wide spectrum of solutions that is absent in 
the case of GR. In particular, there are oscillating solutions [01] i n the case of closed models for 
the Universe, filled by nonlinear scalar field. However, such solutions have nothing to do with 
inflation. 

We conclude therefore that attempts to build nonsingular inflationary cosmological models 
within framework of GR as well as GTG lead to serious complications connected with domination 
of scalar field on the stage of contraction. 

6 Conclusions 

The analysis carried out in this paper shows that the dynamics of the flat cosmological models 
with massive scalar field within the framework of GTG differs significantly from the corresponding 
behavior within GR. In the phase space there are regions when the scalar field acquires imaginary 
part. Corresponding regions are unphysical. The edges of the physical regions are hyperbolas, 
symmetric with respect to the axis if. 

At the stage of expansion most phase trajectories, that start at some point inside the quantum 
edge, i.e. the region, where the energy density is smaller than the planckian value, tend to 
inflationary separatrix. Therefore, inflationary solutions are quite generic feature of cosmological 
models within GTG. 

At the stage of contraction the dynamics of cosmological models within framework of GTG 
changes significantly in comparison with corresponding dynamics within GR due to the presence 
of attractive separatrices. Before reaching the separatrix the effective equation of state of the 
model is p ~ p like in the case of GR. On the separatrix the equation of state tends to the limiting 
case p — p/3. 
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The Hubble parameter at the stage of compression tends to the limiting value H — — 77=, and 
corresponding models have a cosmological singularity in infinite future. Therefore, nonsingular 
inflationary solutions are unstable within the framework of GTG. This conclusion is applicable in 
part to closed and open cosmological models with negligible curvature due to continuity of the 
phase space. 
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